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ABSTRACT 

We calculate the iron line profiles from accretion discs with spiral velocity structures 
around Schwarzschild black holes. We find that quasi-periodic bumps appear in the 
the profiles, thereby providing a test for spiral wave patterns. This study is motivated 
by recent work showing that spiral density waves can result from MHD instabilities 
even in non-self-gravitating discs, and by improved spectral resolution of forthcoming 
X-ray missions. 
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1 INTRODUCTION 

Active galactic nuclei (AGN) and some X-ray binaries are 
thought to be powered by accretion onto a central black hole 
(e.g. Rees 1984). A fluorescent iron Ka emission line with 
rest energy 6.4 keV is observed in the X-ray spectrum of 
many Seyfert galaxies and some X-ray binaries. The line pro- 
file is shaped by Doppler and gravitational redshifts (Fabian 
et al. 1989) and therefore by the geometry and dynamics in 
the innermost regions of the accretion engine (e.g. Fabian et 
al. 2000 for a review). 

The simplest plausible geometry to explain the observed 
profiles is a geometrically thin, optically thick disc with an 
X-ray corona above and below the disc. Some of the X- 
rays emanate directly to the observer, and others are repro- 
cessed by the disc, inducing the iron line. General relativity 
is needed to calculate the redshifts and photon trajectories. 
This has been done for both Schwarzschild (Fabian et al. 
1989) and Kerr (Laor 1991) black holes. The resulting pro- 
files have been successful in reproducing many features of 
the experimental data (e.g. Tanaka et al. 1995; Nandra et 
al. 1997). 

The primary motivation for the present work is that 
MHD instabilities can induce spiral density waves in non- 
self-gravitating accretion discs (Caunt and Tagger 2001, 
Hawley 2001, Tagger and Pellat 1999). It is therefore inter- 
esting to see whether this spiral structure has a significant 
effect on iron line profiles. 

The effect of spiral structure on line profiles has been 
considered before in the context of tidally induced spiral 
shocks in close binary systems (e.g. Chakrabarti and Wiita 
1993, Steeghs and Stehle 1999). The non-axisymmetry of the 
situation was emphasised by Sanbuichi, Fukue and Kojima 



(1994), who also considered special relativistic effects. For a 
review of numerical work on spiral shocks in close binaries 
see Matsuda et al. (2000). 

This work differs from previous considerations of spiral 
structure because the context is iron lines from AGN and 
therefore we use general relativity. For non-self-gravitating 
discs, and discs not in binaries, spiral structure is still possi- 
ble, but requires MHD instabilities (Caunt and Tagger 2001, 
Hawley 2001, Tagger and Pellat 1999). The changes to the 
profile come purely from the spiral velocity structure and 
the corresponding modified redshifts. 

We find that the spiral velocity structure does have a 
significant effect on the profiles. Most notably by introducing 
multiple quasi-periodic sub-peaks into the profile, particu- 
larly visible for tightly packed spirals. For reasonable values 
for the parameters of the spiral structure, the effects on the 
line profiles require better spectral resolution than is cur- 
rently available (e.g. Nandra et al. 1997), but will likely be 
discernable with Astro E-2. We can therefore make predic- 
tions. 

Like warping (see Bachev 1999, Hartnoll and Blackman 
2000) spiral structure induces non-axisymmetric effects on 
the line profiles, though the tighter the winding of the spi- 
ral the less the observed effect. Nevertheless, one immediate 
effect of non-axisymmetry is a time dependence in the line 
profile of rotating discs. Spiral shocks were used in this way 
by Chakrabarti and Wiita (1994) to study variable broad 
optical emission lines in AGN. In fact a time dependence 
in the iron X-ray line is also observed (e.g. Iwasawa et al. 
1996) and has been argued to be problematic for the stan- 
dard disc models (e.g. Weaver and Yaqoob 1998, Sulentic et 
al. 1998). Although the issue remains unresolved, some kind 
of non-axisymmetry of the disc may be playing a role. 
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It is straightforward to adapt the standard calculations 
of line profiles to include spiral structure. The main differ- 
ence is that the azimuthal velocity of the gas is no longer 
Keplerian and the radial velocity is nonzero. Representa- 
tive expressions for the azimuthal and radial velocities are 
given in §2. In §3 we list the formulae needed to compute 
the profiles, in order that our paper be self-contained. In 
particular the redshift formula needs to be generalised for 
non-Keplerian rotation. §4 contains sample profiles and ex- 
planation of the new features. §5 contains a discussion of the 
results. 



2 SPIRAL WAVES IN ACCRETION DISCS 

We use a simple representation of spiral structure in the gas 
velocity, motivated by the numerical results of Caunt and 
Tagger (2001) 
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where ro is a constant parametrising the fall-off with ra- 
dius r, m is the azimuthal wavenumber, k is larger for more 
tightly packed spirals, e is the size of the perturbation, A al- 
lows the velocities to be out of phase. So (e r , e$, m, k, ro, A) 
is the parameter space of discs. We then combine these with 
the angle of disc inclination to the observerm, and the az- 
imuthal viewing angle of the observer. See Fig. 1 for an ex- 
ample disc. Note that in this work we consider spiral waves 
which are stationary over the observation duration. There 
could in addition be a time dependent phase in eqs. (1) and 
(2) to account for a rapid spiral wave propagation. 



3 CALCULATION OF THE LINE PROFILE 
3.1 Total flux 

The total flux observed at frequency v is given by 

F(^) oc J 5{v - u (r))f(r, v)d 2 x, (3) 

where the integral is over the disc, vo(v) is the observed 
frequency of flux emitted from r and /(r, v) is given by 



/(r,^) = 



/ out (r,c(r,n)), 



(4) 



where [^-1 is the relativistic correction to the flux, with 
the frequency of emission, v e — 6.4keV for the flourescent 
iron line. The function / ou t(r, c) is the flux reprocessed and 
re-emitted at r in the direction c, defined such that emit- 
ted photons in this direction from the disc eventually reach 
the observer in the Schwarzschild metric. The asymptotic 
direction to the observer is n. 

The range of emitting region we consider is 6 = ri n < 
r < r out — 60 and < <f> < 2ir, in units with c = G = 
M = 1. The lower bound is the innermost stable orbit of 
the Schwarzschild metric and the upper bound is somewhat 
arbitrary, but unimportant, as the objective is to compare 
our results here with profiles from discs with axisymmet- 
ric density structures that have well known dependences on 
outer radius (e.g. Fabian el al. 1989). Moreover, most of the 



interesting effects will occur in the inner regions. The area 
element has a general relativistic correction, so 

dr 
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(5) 



3.2 Photon paths: direction of emission 

We use paths in the Schwarzschild metric, approximated by 
a first order perturbation. In this approximation, the initial 
direction c of the photon at r reaching the observer, (i.e. in 
asymptotic direction h) is (Hartnoll and Blackman 2001) 



c(r) = -^($ = $ ), 

where $ parametrises the photon path and 
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and the initial value $o is given by 
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C = -2, 
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This perturbative approach is valid throughout most of the 
disc and sufficient for our purposes. Corrections to the pho- 
ton paths due to strong curvature in the innermost regions 
are most significant for very high inclination angle profiles 
(Matt, Perola and Stella 1993). 



3.3 Frequency shifts 

The observed frequency of radiation from r in the 
Schwarzschild metric is (e.g. Hartnoll and Blackman 2001) 



KM 



1 
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(12) 



Here u(r) is the velocity of the emitting material at r, given 
in the previous section. As noted in Hartnoll and Blackman 
(2001), in this context 6 needs to be normalised so that 
£r( 1 r ^ (r) + c (r)c e (r) + c^,(r)c^(r) = jK. Fig. 2 shows the 
redshifts for a sample disc. 



3.4 Reprocessed flux 

We need to calculate / ou t(r,c). We use the standard iron 
line model (Fabian et al. 1989) in which there is assumed to 
be a corona above the central region producing hard X-rays 
which are reprocessed by the disc. The corona is modelled 
by a source which is at height h s directly above the centre 
of the disc. So the flux impinging on the disc at r is 
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(13) 



We will take q = 3 (i.e. 1/r 2 falloff from source) and h e = 10. 
The main objective is to study the effect of spatial inho- 
mogenity on the line profiles, so the exact values of these 
parameters are not important so long as they allow us to 
compare the results with standard line profiles. Assuming 
the re-emission is isotropic, we have 



f out (r,c) oc /i„(r)e« • c, 



(14) 



where e z is the unit normal to the disc. Note that the varying 
density will affect the optical depth but not the total power 
coming from a given radius as long as the optical depth is 
>> 1. We do not incorporate the vertical density structure 
of the disc in calculating our line profiles. 



4 DISCUSSION OF LINE PROFILES 

The region of the parameter space investigated in this sec- 
tion is qualitatively guided by recent numerical results from 
simulations of spiral structure in non-self-gravitating MHD 
accretion discs (Caunt and Tagger 2001, Hawley 2001). 

Fig. 3 shows profiles from discs at 30 degrees inclination 
for various values of k and n>. The spiral velocity structure 
results in the presence of small peaks across the profile. As 
expected, increasing k increases the frequency of the the 
peaks. At lower ro, the spiral structure is more confined 
to the inner regions. Note that e — e r = needs to be 
changed with ro to ensure that the overall normalisation is 
comparable. A concentration of structure at small r reveals 
itself through the prevalence of peaks or bumps appearing 
in the more redshifted or blueshifted regions of the profile. 

Fig. 4 shows profiles from discs at 15 degrees inclina- 
tion for various k and a profile with no spiral structure for 
comparison. A larger range of k is considered in comparison 
to Fig. 3. Again the correlation between number of peaks 
and k is evident. 

Fig. 5 shows profiles from discs at an inclination angle 
of 60 degrees. The profiles from discs at high inclination 
angles are naturally more spread out. The effects of k and 
ro are as before. It is particularly noticeable in these profiles 
that the sub-peaks occur only outside of the standard double 
peaks. This is due to the falloff with r; the spirals are more 
prominent in the inner regions of the disc, which are more 
Doppler shifted. The atypical blueshifted region of the k = 
1.2, ro = 15 graph is also interesting, the sub-peaks evolve 
into a step-like form. 

All discs considered above have m = 1, A = 0, e r — e$ 
fixed given ro, and the observer is at the same azimuthal 
angle in all cases. Varying e with ro fixed has the obvious ef- 
fect of increasing the relative magnitude of the perturbations 
to the usual line profile from discs without spiral structure. 
Varying m, A or the azimuthal angle of the observer does 
not appear to have a significant qualitative effect on the 
profiles. 



red or blue-shifted regions, indicating spiral structure in the 
inner disc. The profiles can be compared with data to infer a 
constraint on dynamical models predicting spiral structure. 
Figs. 3, 4 and 5 show what sort of values of e, for a given ro, 
are consistent with present data. 

Though not easily constrained with current instru- 
ments, this quasi-periodic bumpiness should be resolvable 
with Astro E-2, which carries the X-ray Spectrometer (XRS) 
of spectral resolution AE/E ~ 0.002 in the 6.4keV range 
(Kelley et al. 1999) which should be sufficient for a wide 
range of k. For very low k, there are fewer sub-peaks, and 
their exact shape and position depend on the variables 
m, A, ^- and the azimuthal angle of the observer. 
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5 CONCLUSION 

The most interesting overall feature in our line profiles is the 
presence of multiple, quasi-periodic sub-peaks in the most 
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Figure 4: Clockwise from top left: The first three are profiles 
from discs at inclination angle f5 degrees with m = 1, A = 0, 
r = 15,e = 0.2. From left to right: k = 0.4, k = 1.2, k = 2.0. The 
profile at the bottom left is from a disc with no spiral structure, 
for comparison. 
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Figure 5: As for Fig. 3 but at inclination angle of 60 degrees. Left 
column: k = 0.4. Right column: k = 1.2. Top row: ro = 5, e = 0.3. 
Bottom row: ro = 15, e = 0.2. 

Figure 1: The spiral structure of the radial velocity in a 
disc with m — 1, k — 1.2, ro = 15. 

Figure 2: The redshifts of the disc of Fig. 1, viewed at 
an inclination angle of 30 degrees. Darker shading indicates 
higher redshift. The disc is rotating anticlockwise, the ob- 
server is situated to the right. 

Figure 3: Profiles from discs at inclination angle 30 degrees 
with m = 1, A = 0. Left column: k — 0.4. Right column: 
k — 1.2. Top row: ro = 5, e = 0.3. Bottom row: ro = 15, e = 
0.2. 
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